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We report the results of a search for D -D mixing in D° — > K~*~7i~ decays based on 400 ftr 1 of 
data accumulated by the Belle detector at KEKB. Both assuming CP conservation and allowing for 
CP violation, we fit the decay-time distribution for the mixing parameters x' and y' , as well as for 
the parameter Rd, the ratio of doubly-Cabibbo-suppressed decays to Cabibbo-favored decays. The 
95% confidence level region in the (x' 2 ,y') plane is obtained using a frequentist method. Assuming 
CP conservation, we find x' 2 < 0.72 x 10" 3 and -9.9 x 10~ 3 < y < 6.8 x 10~ 3 at the 95% confidence 
level; these are the most stringent constraints on the mixing parameters to date. The no- mixing point 
(0, 0) has a confidence level of 3.9%. Assuming no mixing, we measure Rd = (0.377±0.008±0.005)%. 

PACS numbers: 13.25.Ft, 11.30.Er, 12.15.Ff 



The phenomenon of mixing has been observed in the 
K°-~K° and B°-B° systems, but not yet in the D a -D a 
system. The parameters used to characterize mixing are 
x = Am/r and y = Ar/(2T), where Am and AT arc 
the differences in mass and decay width between the two 
neutral D mass eigenstates, and T is the average width. 
The mixing rate within the Standard Model is expected 
to be small 0. The largest predicted values, including 
long-distance effects, arc of order |x| < \y\ ~ (1CP 3 — 
10~ 2 ), and are reachable with the current experimental 
sensitivity. Observation of |ar| 3> \y\ or CP violation 
(CPV) in D°-D° mixing would constitute unambiguous 
evidence for new physics. 

The "wrong-sign" (WS) process, D° — > K + ir~, 
can proceed cither through direct doubly-Cabibbo- 
suppressed (DCS) decay or through mixing followed by 
the "right-sign" (RS) Cabibbo-favored (CF) decay D° -> 
The two decays can be distinguished by 
the decay-time distribution. For |a;|,|y| <C 1, and as- 
suming negligible CPV, the decay-time distribution for 
D° — > K + it~ can be expressed as 
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~dt 



oc e 



-rt 



(1) 



where Rd is the ratio of DCS to CF decay rates, x' = 
x cos 5 + y sin 5, y' = y cos 5 — x sin 5, and S is the strong 
phase difference between the DCS and CF amplitudes. 
The first (last) term in brackets is due to the DCS (CF) 



amplitude, and the middle term is due to interference be- 
tween the two processes. The time-integrated rate (Rws) 
for D° -> K+tt- relative to that for D° -> R-tt + is 

R D + VR^y' + {x ,2 + y' 2 )/2. 

To allow for CPV, we apply Eq. © to D° and D° sep- 
arately. This results in six observablcs: x' +2 , y' + } 
for D° and {R}j,x' 2 ,y' } for D° . CPV is parameter- 

R- D )/{R + D + R D ) 
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ized by the asymmetries Arj 
and A M = {Rjj - R M 

(x' ±2 +y /±2 )/2. A D and A M characterize CPV in DCS 
decays and in mixing, respectively. The observables are 
related to x' and y' via 
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(y cos <j> =F x sin 4>) , (3) 



where is a weak phase and characterizes CPV occur- 
ring in interference between mixed and unmixed decay 
amplitudes. To avoid ambiguity, we restrict <j) to the 
range \<j>\ < ir/2. 

This method has been exploited in previous studies 
Ellil 00 in our previous measurement based on a 90 
fb _1 data sample, the value of y' was found to be slightly 
positive although compatible with zero 0- Here we ex- 
ploit the much larger data set now available to search for 
D°-D° mixing with significantly higher sensitivity. 
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In this Letter we present improved results of an anal- 
ysis of 400 fb _1 of data, setting more stringent lim- 
its on mixing and CPV parameters. The data were 
recorded by the Belle detector at the KEKB asymmetric- 
energy e + e~ collider The Belle detector |tj includes 
a silicon vertex detector (SVD), a central drift chamber 
(CDC), an array of aerogel threshold Cherenkov counters 
(ACC), a barrel-like arrangement of time-of- flight scintil- 
lation counters (TOF), and an electromagnetic calorime- 
ter. The first 157 fb _1 of data were taken with a 2.0 cm 
radius bcampipe and a 3-layer SVD, while the subsequent 
243 fb _1 were collected with a 1.5 cm radius beampipe, 
a 4-layer SVD, and a small-cell inner drift chamber |10j . 

We reconstruct D° candidates from the decay chain 
D*+ -> 7T+£>°, D° -> K ± nT. Here, tt s denotes the low- 
momentum (slow) pion, the charge of which tags the fla- 
vor of the neutral D at production. We select D° can- 
didates by requiring two oppositely-charged tracks, each 
with at least two SVD hits in both r-<f> and z coordinate, 
satisfying K and 7r identification selection criteria. These 
criteria are C > 0.5 for K and C < 0.9 for n, where C is 
the relative likelihood for a track to be a K based on the 
response of the ACC and measurements from the CDC 
and TOF. These criteria have efficiencies of 90% and 
94%, and ir/K misidentification rates of 10% and 17%, 
respectively. To reject background candidates in which 
the K is misidentified as n and the ir is misidentified as 
K, we recalculate tuk-k with the K and ir assignments 



(swapped) 



< 28 



swapped and reject events with 

MeV/c 2 (~ 4.5c). A D* + candidate is reconstructed by 
combining a D° candidate with a ir s candidate; the re- 
sulting D* + momentum in the e + e _ center-of-mass frame 
{pd*) is required to be > 2.7 GeV/c in order to eliminate 
BB events and suppress the combinatorial background. 

The D° vertex is obtained by fitting its daughter 
tracks. The D* vertex is taken as the intersection of the 
D trajectory with the interaction region. We constrain 
7r s to originate from the obtained D* vertex. A good 
X 2 for each vertex fit is required. The D° proper decay 
time t is then calculated. We require the uncertainty 
of the decay-time a t to be less than 0.7 ps (typically, 
(74^0.13 ps). 

The selection criteria for particle identification, the 
X 2 of vertex fits, and pn* are obtained by maximizing 
N S i g /y/N s i g + A^bkg, where iV s i g (iVbkg) is the expected 
number of WS signal (background) events estimated from 
data in the RS signal (WS sideband) region. We assume 
^?ws = 0.37% in the calculation. The optimized val- 
ues of the selection criteria are similar to those used pre- 
viously 0. 

We select events satisfying 1.81 GeV/c 2 < rriKir < 1-91 
GeV/c 2 and < Q < 20 MeV, where Q = m Kn „ s - 
rriKir — is the kinetic energy released in the decay. 
About 5% of selected events have two or more D* can- 
didates associated with a single D° candidate. If these 
D* candidates have opposite sign, the event is rejected; 



this reduces random n s background (see below) by 30% 
while reducing the signal by only 1%. If the D* candi- 
dates have the same sign, then we choose the candidate 
that has the best x 2 resulting from the vertex fit. 

We determine RS and WS event yields from a two- 
dimensional fit to the rriK-n-Q distribution. There are 
four significant background sources in the WS sample: 
(a) random n s background, in which a random tt + is 
combined with a U° — > K + n~ decay; (b) D* + — > D°7r+ 
followed by D° decaying to > 3-body final states; (c) 
decays; and (d) combinatorial. They are denoted 
as rnd, d3b, ds3 and cmb in turn. These background 
shapes are obtained from Monte Carlo (MC) simulation 
and fixed in the fit. When fitting the RS sample, the 
parameters for the signal shape are floated; when fitting 
the WS sample, these parameters are fixed to the values 
obtained from the RS fit. We find 1073993±1108 RS and 
4024 ± 88 WS signal events, and the ratio of WS to RS 
events is (0.375±0.008)% (statistical error only). The ra- 
tio of WS signal to background is 1.1, about 20% higher 
than that of our previous study 0. The background is 
composed mostly of random 7r s (51%) and combinatorial 
(35%) events. Figure shows the rriK-n and Q distribu- 
tions superimposed with projections of the fit result. The 
WS projections for ttikt: and Q are shown for a 3a signal 
interval in Q and mx^, respectively. The contribution of 
decays is too small to be seen. 
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FIG. 1: The distribution for (a) RS m K -n with < Q < 20 
MeV; (b) RS Q with 1.81 GeV/c 2 < m K7T < 1.91 GeV/c 2 ; (c) 
WS mx, with 5.3 MeV< Q < 6.5 MeV; and (d) WS Q with 
1.845 GeV/c 2 < m & < 1.885 GeV/c 2 . Superimposed on the 
data (points with error bars) are projections of the tuktv-Q 
fit. 

The decay-time fitting procedure is similar to that of 
our previous measurement 0] but with several improve- 
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ments as discussed below. We determine Rd, x' 2 and 
y' by applying an unbinned maximum likelihood fit to 
the distribution of WS proper decay time, considering 
the 4<r region \m K ^ — m D a\ < 22 MeV/c 2 and \Q — 5.9 
MeV| < 1.5 MeV. We determine background shapes by 
fitting events in an mjf, sideband (this contains no signal 
or random 7r s events). The probability density function 
(PDF) for the WS signal is given by Eq. denoted as 
P s ig, convolved with a resolution function R s i g (t). The 
latter is represented by a sum of three Gaussians with 
widths <jj = Sj ■ a t (J = 1 — 3) and a common mean 
(the error a t varies event-by-event). The decay-time dis- 
tributions for the backgrounds from random n s , 3-body 
and Dt-, decays are exponential, Pk = e~'/ Tfc (£;=rnd, 
d3b, ds3), while the distribution of combinatorial back- 
ground is taken to be a Dirac delta function 5(t). The 
distributions are also convolved with the corresponding 
resolution functions Rk (fc=d3b, ds3, cmb) which depend 
on a t . For the main background contribution of random 
ir s , the resolution function and the lifetime are the same 
as those of the signal, since the ir s does not affect the D° 
vertex reconstruction. We define a likelihood value for 
each (ith) event as a function of Rd, x' 2 and y': 

poo 

Pi= / dt> [{r sig P sig (t';R D ,x' 2 ,y') 
Jo 

+ /rnd^rnd(i') } Rsig(U - t ) + /d3b-Pd3b (,t')Rd3b (*i — t') 

+ fds 3 Pds3(t')R ds3 (U - f) +f> mh 6(t')R cmh (U - f )] . 

(4) 

Here, the fractions =sig, rnd, d3b, ds3 or cmb) are 
determined on an event-by-event basis as functions of 
m&, Q and a t . 

The fitting procedure is implemented in steps as fol- 
lows. First we fit the RS sideband region using a simple 
background model to obtain parameters of i?d3b- Then 
we fit the same events using a full background model 
as in Eq. Q, which yields i? C mb and Td3b for RS back- 
ground. We fit the RS signal region with these back- 
ground parameters fixed, and obtain parameters of R slg 
(the scaling factors Sj, and the mean value and frac- 
tions of the individual Gaussians) and the D° lifetime 
r D a- The latter is found to be 409.9 ± 0.7 fs, in good 
agreement with the world average value [ll|. The % 2 of 
the fit projection on the decay-time distribution is 64.0 
for 60 bins. We use different resolution parameters for 
the two SVD configurations. We then fit the WS sam- 
ple. We fit the WS sideband region with i?d3b fixed from 
the RS sideband fit and the Df \ contribution fixed from 
MC calculations; this yields i? cm b and Td3b for WS back- 
ground. Finally, we fit to the WS signal region with these 
background parameters, i? s ig and t d o fixed. In the final 
fit, Rd, x' 2 and y 1 are the only free parameters and are 
determined by maximizing the extended log-likelihood 
function ln£ = X^ m ^ + m £-R. The function Cr is 
a Gaussian that constrains the ratio Rws(Rd, x' 2 , y') to 



be near the value obtained from the ttlktv-Q fit; this is 
needed because P s ig(t'; RD,x l2 ,y') is normalized to unity. 

The main improvements in the decay-time fitting pro- 
cedure with respect to that of our previous measurements 
on a smaller data set consist of using an improved res- 
olution function and optimized coefficients F° r the 
latter, we include a dependence on at, as this variable 
substantially improves the discrimination between signal 
decays and combinatorial background. We determine the 
at distribution for combinatorial background by fitting 
WS data. To check the correctness of this method we 
generate MC samples with the same size as data, add 
the corresponding amount of backgrounds, and repeat 
the fitting procedure. For a wide range of (x' 2 ,y') values, 
the fit recovers the input values well within the statistical 
uncertainty. If the at dependence is not included in 
the fit obtains values shifted by 0.5-1 statistical standard 
deviation with respect to the input values. 



TABLE I: Summary of results including systematic errors. 



Fit case 


Parameter 


Fit result 95% C.L. interval 
(xlO -3 ) (xlO -3 ) 


No CPV 


Rd 

y' 

Rm 


3.64 ±0.17 (3.3, 4.0) 
0.18ig;g < 0.72 

0.6+39 (- 9 - 9 . 6 - 8 ) 

(0.63 x 10" 5 , 0.40) 


CPV 


A d 


23 ±47 (-76, 107) 




A m 


670 ± 1200 (-995, 1000) 




x' 2 


< 0.72 




>/ 


(-28, 21) 




Rm 


< 0.40 


No mixing 


Rd 


3.77 ± 0.08(stat.) ± 0.05(syst.) 




-2000 2000 4000 

Proper time (fs) 



FIG. 2: The decay-time distribution for WS events satisfying 
\m Kw - m D o \ < 22 MeV/c 2 and \Q - 5.9| < 1.5 MeV. Super- 
imposed on the data (points with error bars) are projections 
of the decay-time fit when no CPV is assumed. The mixing 
and interference terms are shown at the 95% confidence level 
upper limit (95% UL) for mixing. 

Tablc[I]lists the results from three separate fits. For the 
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first fit, we require CP to be conserved. The projection 
of this fit superimposed on the data is shown in Fig. |5J 
the x 2 of the projection is 54.6 for 60 bins. The central 
value of x' 2 is in the physically-allowed region x' 2 > 0. 
The correlation between x 12 and y' is —0.909. The results 
for the two SVD subsamples are consistent within 0.6 a. 
For the second fit, we allow CPV and fit the WS D° and 
D° samples separately. We calculate Ad and Am (see 
Table Pi, and solve for x' 2 , y' and <f> using Eqs. J2J) and 
©. We obtain \<j>\ = (9.4 ± 25.3]° or (84.5 ± 25.3)° for 
the same or opposite signs of x' and x' . Finally, for 
the last fit we assume no mixing and set x' 2 = y' — 0. 

We apply the method described in Ref. to obtain 
the 95% confidence level (C.L.) region and take into ac- 
count the systematic errors. Figurc|21shows the 95% C.L. 
contours with and without CPV allowed. For the case 
of no CPV, the allowed area of (x' 2 ,y') values is smaller 
than that of our previous measurement by a factor of 2.2. 
The CPV contour has a complicated shape due to there 
being two solutions for (x', y') when solving Eqs. and 
0, depending on the signs of x ,:lz . 
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FIG. 3: 95% C.L. regions for (a;' 2 , y'). The point is the best fit 
result assuming CP conservation. The dotted (dashed) line 
is the statistical (statistical and systematic) contour for no 
CPV . The solid line is the statistical and systematic contour 
in the CPV- allowed case. 

We evaluate systematic errors by varying parameters 
used to select and fit the data within their uncertainties. 
The sources of systematic error include event yields and 
imperfect modeling of backgrounds and uncertainties in 
the decay-time PDF's. The former were estimated by 
changing the selection criteria (kaon and pion identifica- 
tion, x 2 °f vertex fits, and the D* momentum) and thus 



the signal to background ratio over a significant range. 
The significance of an individual systematic shift is found 
by calculating m 2 = — 2[lnL(a now ) — lni(ao)]/2.3, with 
<^ncw — (^ new; 2/'new) denoting the result of the fit with 
the modified parameter and <So the result from the de- 
fault fit. The factor 2.3 corresponds to 68% confidence 
in two dimensions. The largest shift occurs for the D* 
momentum selection; it is found to be m 2 = 0.083. The 
parameters of functions fitted to the vtik-k and Q distribu- 
tions were also varied by their corresponding uncertain- 
ties and the decay-time fit was repeated. The resulting 
systematic error is found to be small. The influence of 
a t on the fractions fl is checked by obtaining the com- 
binatorial background at PDF from the fit to sideband 
events. Repeating the time fit with the modified fj. yields 
m 2 = 0.030. The same value is found when varying all of 
the fixed parameters entering the decay-time PDF's by 
their uncertainty. Adding in quadrature the significances 
of all shifts due to possible systematic uncertainties, we 
find the overall scaling factor yl + X] m 2 = 1-12. We 
increase the 95% C.L. statistical contour by this factor 
to include systematic errors. 

We show the contour with systematic errors included 
in Fig.OHas a dashed line in the CP-conserving case and 
as a solid line in the general case. In the case of no CPV , 
the no-mixing point x' 2 — y' — lies just outside the 95% 
C.L. contour; this point corresponds to 3.9% C.L. with 
systematic uncertainty included. The two-dimensional 
95% C.L. intervals of parameters listed in Tablc|Uare ob- 
tained by projecting these contours onto the correspond- 
ing coordinate axes. In the case of CPV, because the 
95% C.L. contour includes the point x' 2 = y' = 0, we 
cannot constrain <f> at this confidence level. 

In summary, we have searched for D°-D° mixing and 
CP violation in "wrong-sign" D° — * K + ir~ decays using 
a 400 fb _1 data sample. Assuming negligible CP viola- 
tion in the D° system, we obtain x' 2 < 0.72 x 10~ 3 and 
-9.9 x 10~ 3 < y' < 6.8 x 10~ 3 at 95% C.L. These re- 
sults supercede our previous measurement and represent 
the most stringent limits on D°-D° mixing parameters 
to date. The data exhibits a small preference for positive 
x' 2 and y'\ the no- mixing point x' 2 = y' = corresponds 
to a C.L. of 3.9%. 
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